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Abstract 

The two body charmless decays of B s meson to light vector mesons are analyzed within 
the framework of QCD factorization. This approach implies that the nonfactorizable 
corrections to different helicity amplitudes are not the same. The effective parameters 
for helicity h = 0, +, — states receive different nonfactorizable contributions and hence are 
helicity dependent, contrary to naive factorization approach where are universal and 
polarization independent. The branching ratios for B s — > VV decays are calculated and 
we find that branching ratios of some channels are of order 10 -5 , which are measurable at 
future experiments. The transverse to total decay rate IY/r is also evaluated and found 
to be very small for most decay modes, so, in charmless B s — ► VV decays, both light 
vector mesons tend to have zero helicity. 
PACS Numbers 13.25.Hw, 12.38. Bx, 12.15Mm 
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1 Introduction 



The charmless two-body B decays play a crucial role in determining the flavor parameters, 
especially the Cabibbo-Kobayashi-Maskawa(CKM) angels a, (3 and 7. With precise measure- 
ments of these parameters, we can explore CP violation which is described by the phase of 
the CKM matrix in the standard model(SM). Recently there have been remarkable progresses 
in the study of exclusive charmless B decays, both experimentally and theoretically. On the 
experimental aspect, many two-body non-leptonic charmless B decays have been observed by 
CLEO and S-factories at KEK and SLAC Ej and more B decay channels will be measured 
with great precision in the near future. With the accumulation of data, SM can be tested in 
more detail. Theoretically, several novel methods have also been proposed to study the non- 
factorizable effects in the hadronic matrix elements, such as QCD factorization(QCDF) jH], the 
perturbative QCD method(pQCD) |lj and so on. Intensive investigations on hadronic charmless 
two-body B u d decays with these methods have been studied in detail 13 El E] • 

The extension of QCDF from B u ^ decays to B s decays has also been carried out by several 
authors ^01^3- m principle, the physics of the B s two-body hadronic decays is very similar to 
that for the B^ meson, except that the spectator d quark is replaced by the s quark. However, 
the problem is that B s meson oscillates at a high frequency, and nonleptonic B s decays have 
still remained elusive from observation. Unlike the B u ^ mesons, the heavier B s meson cannot 
be studied at the B-factories operating at the T(4s) resonance. But it is believed that in the 
forthcoming hadron colliders such as the Collider Detector at Fermilab (CDF), DO, DESY ep 
collider HERA-B, BTeV, and CERN Large Hadron Collider (LHCb), CP violation in B s system 
can be observed with high accuracy. This makes the search for CP violation in the B s system 
decays very interesting. 

In the papers [T2|, the authors have studied systematically the B s — > PP,PV de- 
cays(here P, V denote pseudoscalar and vector mesons respectively) with QCD factorization, 
and intensive phenomenological analysis has been made. Since the B — > VV modes reveal 
dynamics of exclusive B meson decays more than the the B — > PP and PV modes through the 
measurement of the magnitudes and the phases of various helicity amplitudes, in the present 
work we plan to make a detail study of B s — > VV decays within the same framework of QCD 
factorization. We find that, contrary to the generalized factorization approach |10j. nonfactor- 
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izable corrections to each helicity amplitude are not the same; the effective parameters vary 
for different helicity amplitudes. The transverse to total decay rate IV /T is very small for most 
decay modes, so in the heavy quark limit, both light vector mesons in charmless B s — > VV 
decays tend to have zero helicity. Branching ratios for some decay modes are found of order 
10~ 5 , which could be measured at LHCb. 

This paper is organized as follows. In Sec. II, we outline the necessary ingredients of the 
QCD factorization approach for describing the B s — ► VV decays and calculate the effective 
parameters af. Input parameters, numerical calculations and results are presented in Sec. 
III. Finally we conclude with a summary in Sec. IV. The amplitudes for charmless two-body 
B s — > VV decays are given in Appendix. 

2 B s — > VV in QCD factorization approach 
2.1 The effective Hamiltonian 

Using the operator product expansion and renormalization group equation, the low energy 
effective Hamilization relevant to hadronic charmless B decays can be written as ^3] 

H eff = ^[A^dOr + C 2 O u 2 ) + A c (dOJ + C 2 O c 2 ) 

10 

-x t CE + c 7l o 7l + c 8g o 89 )} + h.c (i) 

i=3 

where Aj = VibV* q are CKM factors and Cj(/x) are the effective Wilson coefficients which have 
been reliably evaluated to the next-to-leading logarithmic order. The effective operators Oi can 
be expressed as follows: 



0? 


= (ub) v _ A (qu)v_ A , 


O u 2 


= (Uab^V-AtqpUajv-A, 


0{ 


= (cb) v _ A (qc) v _ A , 


oi 


= (c a b/ 3 ) v _ A (q f3 c a )v-A, 


°3(5) 


= {qb)v-AH (q'q') v-A(v+A) , 
q' 


^4(6) 


= (qabp)v-AY,(q'(3q' a )v-A(V+A), 

q' 


^7(9) 


= I(qb)v-Aj2e q '(q'q') v +A(y-A), 
q' 


Os(io) 


= I(qabp)v-Aj2e q '{q' p q' a )v+A(v-A) 
q' 




= g^m^cr^l + 7 5 )&a^, 




= ^m b q a a^(l + l5 )T^bpG* u . 


Where q 


= d,s and q' denotes all the active quarks 


at the scale \x = 0(mb), i.e., q' = u, d, 



(2) 
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2.2 The factorizable amplitude for B s — > VV 

To calculate the decay rate and branching ratios for B s — > VV decays, we need the hadronic 
matrix element for the local four fermion operators 

(V^V^iM^v-Aiq^v-AiBs), (3) 

where Ai, A2 are the helicities of the final-state vector mesons V\ and V2 with four-momentum 
Pi and P2, respectively. In the rest frame of B s system, since B s meson has spin zero, we have 
Ai = A2 = A. Let X( BaVl ' V2 ' denote the factorizable amplitude with the vector meson V2 being 
factored out, under the naive factorization(NF) approach, we can express as 

x (B aVl ,v 2) = (y 2 |(g 2g3 ) y _ A |o}(T/ 1 |(g 1 6)^|5 s ). (4) 
In term of the decay constant and form factors defined by EE 

(V(p,e*)\qi,q'\0) = -if v m V £^ (5) 

< V(p,e*)\qi»(l - ls)b\B s {p B ) > = -e;(m B + m v )A B ° v (q 2 ) + (p B +p),(e* ■ p B ) }^ 

m B + my 

+q,(e*-PB) 2 -^[A^ v ^)-Ar(q 2 )] 

-ie^e^ (6) 
m B + my 

where q = p B — p and the form factors obey the following exact relations 

A 3 (0) = A (0), 

Af-Vtf) = mB + mv A B ° v (q 2 ) - mB ~ mv A B ° v (q 2 ). (7) 
2mv 2uiv 



With above equations, the factorizable amplitude for B s —>■ V1V2 can be written as 

— T 2A BsVl (m 2 ) 

x (b sVi ,v 2 ) = . f (el-£* 2 )(m B .+m v Mf' V KK l )-(el-pBM-PB) 2 j Va) 

m Ba + m Vl 

(8) 



+ieiwa/3£2 £i PbP ; 

m Ba + m Vl _ 

where p B (m Bs ) is the four-momentum(mass) of the B s meson, m^^) and m^^) are the 
masses(polarization vectors) of the two vector mesons V\ and V2 respectively. Here and in the 
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following throughout the paper we use the sign convention e 0123 = — 1. Assuming the V1CV2) 
meson flying in the plus(minus) z-direction carrying the momentum pi(p2), we get 

-i/v. 



X (B s Vi,V 2 ) 



2my 1 



m Ba +my 1 



m B s ~ m \ ~ m v 2 )( m B s + m Vl )Af aVl (m 2 V2 ) 



l 2 *»'V 2 > 



-if V2 m V2 [{m Ba + m Vl )Af sVl (K 2 ) ± ^p^V B ^{p^)\ ee h ± for A = ±, 



(9) 

where A = 0, ± is the helicity of the vector meson and p c = \pi\ = \p 2 \ is the momentum of 
either of the two outgoing vector mesons in the B s rest frame. 

In general, the B s — > V1V2 amplitude can be decomposed into three independent helicity 
amplitudes H , H + and H_, corresponding to A = 0, + and — respectively. We use the notation 

H\=< Vi(X)V 2 (X)\HefffB s > (10) 

for the helicity matrix element and it can be expressed by three independent Lorentz scalars 
a, b and c. The relations between them can be written as [HI E| 

Hx = e%et (ogT + —^—PbPb + —^e^p la p 2 ^ , (11) 
F \ m Vl m V2 m Vl m V2 J 

where the coefficient c corresponds to the p-wave amplitude, a and b to the mixture of s- and 
d-wave amplitudes. The helicity amplitudes can be reconstructed as 

Hq = -ax -b(x 2 -1), (12) 



H± = -aTcVx 2 - 1, (13) 

where x = j^-^pj- Given the helicity amplitudes, the decay rate and the branching ratio for 
B s — > V1V2 can be written as 



FiBs^V^) = -^(\H \ 2 +\H + \ 2 +\H. 
87rm z B v 



|2 



S, 



Br{B s ^V x V 2 ) = t Bs ^^(\H \ 2 +\H + \ 2 + \H„\ 2 )s, (14) 



with s — 1/2 for two identical final states and s = 1 for the other cases, where r Bs is the 
lifetime of the B s meson, and p c is given by 



2.3 QCD factorization for B s VV 

Under the naive factorization(NF) approach, the coefficients a, are given by dj = Cj + ^Cj + i 
for odd i and aj = Cj + ^r-Cj_i for even i, which are obviously independent of the helicity A. In 
the present paper, we will compute the nonfactorizable corrections to the effective parameters 
a% , which however are not the same for different helicity amplitudes Ho and H±. 

The QCD-improved factorization(QCDF) approach advocated by Beneke et al. [3] allows 
us to compute the nonfactorizable corrections to the hadronic matrix elements (ViV2\Oi\B s ) in 
the heavy quark limit, since in the — > oo limit only hard interactions between the (B s Vi) 
system and V2 survive. In this method, the light-cone distribution amplitudes(LCDAs) play an 
essential role. Since we are only concerned with two light vector mesons in the final states, the 
LCDAs of the light vector meson of interest in momentum configuration are given by fSJ E] 

Mt = Mj all +Ml ± (16) 

with (here we suppose the vector meson moving in the n_-direction) 



Ml = A i- •!(«) - ^ 



l(a)V / \ 

gy {u) + ie^ pa e ± nZ.n + 'f'y s 



(17) 

where n± = (1, 0, 0, ±1) are the light-cone null vectors, u is the light-cone momentum fraction 
of the quark in the vector meson, fy and fy are vector and tensor decay constants, and E is the 
energy of the vector meson in the B s rest system. In Eq . ()17jl . <f>Y (u) and are leading-twist 

distribution amplitudes(DAs), while g±^ V (u) and g'^f 1 ^ (u) = dg± du ^ are twist-3 ones. Since 
the twist-2 DA contributions to the vertex corrections and hard spectator interactions 

vanish in the chiral limit, and furthermore, the contributions of the twist-3 DAs ^'^(u) are 
power suppressed compared to that of the leading twist ones for the helicity zero case, therefore 
we will work to the leading-twist approximation for longitudinally polarized states and to the 
twist-3 level for transversely polarized ones. We note that the same observation has been made 
by Cheng and Yang [TH] in studying B u d — > (f>K*. 

In the heavy quark limit, the light-cone projector for B meson can be expressed as |5j ITH] 

M% = -^p- [(1+ ^) 75 {$f(0+ (18) 
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with the normalization condition 



Jo Jo 



(19) 



where £ is the momentum fraction of the spectator quark in the B meson. 

Equipped with these preliminaries, we can now calculate the nonfactorizable corrections to 
the effective parameters a% systematically. After direct calculations, we get 



no 



C 2 + 

c 3 + 
c 4 + 



c 2 


a s C F 


N c 


AttN c 


Ci 


a s C F 


N c 


AttN c 


c A 


a s C F 


N c 


AttN c 


c 3 

N c 


a s C F 
AttN c 



AttN, 



c 



CaCtf + /&), 
Ci (fi+fn), 

cm + m 



(C 3 --C 9 )[G h (s q ) + G h (s b ) 



vl 



A 
A, 



+ (C 4 + C 6 )J2G h (s l ) + -(C 8 



C 10 )^G h { Si ) + C 8g Gl 



C 5 + 

c 6 + 
c 7 + 



9l 

N c 

Nc 

N c 
CV 

Nc 



Ot-s C F 
Att~W 



c 



Ot-s C F 



c 



N 



c 



AttN c 



C 
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C 9 a s C F h h 
7Vc + 4^7Vc C9(// +///) ' 



(20) 



where C F = (N c — l)/2Nc, and Nc = 3 is the number of colors, Si = mf/ml and q = 
d, s (determined by the b — > d or b — > s transition process). The superscript h denotes the 
polarization of the vector meson(which is equivalent to A, but for convenience we shall adopt 
h in the following) where h = denotes helicity state and h = ± for helicity ± ones. In the 
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expression a\, the upper value in parenthesis corresponds to h — state, while the lower value 
to h = ± ones. 

In Eg. (120)1 . fj denotes the contributions from the vertex corrections. In the naive dimen- 
sional regularization(NDR) scheme for 75, it is given by 



/? = -121og-^-18+ / du$Y 2 (u)(3- — — logu-3i7r), 

Jo 11 1 — it 

= -mog^-ie+fdrf^T^^^ 

m& jo 4 L 1 — it 



+2 / dxdy[ y - — =F \ > , 21 

jo xy xu + y y{xu + y) J 

where ( = +1 or —1, corresponding to (V — A) ® (V — A) or (V — A) <S> (V + A) current 
respectively. It is obvious that ff has the same expression as the hard scattering kernel Fm 2 
for B — > nix mode [HUE] as it should be. 

For hard spectator interactions, supposing V\ be the recoiled meson and V2 the emitted 
meson, explicit calculations for fjj yields 

ff, = jps^ u*m rJM r d /i"w 



Afc ft, jo £ Jo v Jo u 

± = £ttJto 2(1 ± 1) / d? £M / d „?M / du(9 <* (u) T £& 

jVc m Bs h± Jo 4 jo ?r jo 4 



4 U m; 2 



(rH^o^, (22) 



with v — 1 — v, and ho,h± given by Eq. Q. In Eq. (|2~2j) . when we adopt the asymptotical 
form for the vector meson LCDAs, there will be a logarithmic infrared divergence with regard 
to the v integral in ffj, which implies that the spectator interaction is dominated by soft gluon 
exchanges in the final states. In analogy with the treatment in works [HJ El IHj , we parameterize 
it as 

X h = f 1 dx- = \og^(l + p H e^), (23) 
Jo x Ah 

with (ph, <Ph) related to the contributions from hard spectator scattering. Since the param- 
eters (ph,4>h) are unknown, how to treat them is a major theoretical uncertainty in the 
QCD factorization approach. In the later numerical analysis, we shall take = 0.5GeV, 
(ph, (f>h) = (0, 0) [TT] as our default values. 



In calculating the contributions of the QCD penguin-type diagrams, we should pay attention 
to the fact that there are two distinctly different contractions argued in jH] . With this in mind, 
the nonfactorizable corrections induced by local four-quark operators Oi can be described by 
the function G h (s) which is given by 

G°(s) = -? + !logA_ 4 f'du^f^g^s), 
6 6 mh Jo 11 

G^s) = -| + |log^-2/ 1 dn(^H T ^^Mn, S ), (24) 
6 6 mh Jo 4 



with the function 



g(u,s)— / da; xxlog [s — xx(l — u) — ie]. (25) 
Jo 



In Eq. fl20|) . we also take into account the contributions of the dipole operator Og 5 which 
will give a tree-level contribution described by the function G 1 ! defined as 



G° = 


,i 2$j*(u) 
du 

Jo 1 — u 




C 9 + = 


Jo 


4 


G 9 = 


['du (g { ^(u)- 
Jo 


f(a)V 2 1 \ 

4 



u 



(26) 



Due to < V|gig 2 |0 >= 0, B s — > V1V2 decays do not receive nonfactorizable contributions 
from ag and ag penguin terms as shown in Eq. (j20jh 



3 Numerical results and discussions 

To proceed, we use the next-to-leading order Wilson coefficients in the NDR scheme for 75 

Ci = 1.078, C 2 = -0.176, C 3 = 0.014, C A = -0.034, C 5 = 0.008, C e = -0.039, 
C 7 /a = -0.011, C 8 /a = 0.055, C 9 /a = -1.341, C 10 /a = 0.264, C 8g = -0.146. (27) 

at /i = rrih = 4.66 GeV, with a being the electromagnetic fine-structure coupling constant. 
For quark masses, which appears in the penguin loop corrections with regard to the functions 
G h (s), we take 

m u = md = ms = o, m c = 1.47 GeV, m b = 4.66 GeV. (28) 
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As for the CKM matrix elements, we adopt the Wolfenstein parametrization up to 0(X 3 ) 

( ~\ r T/ T/ \ / 1 \2/o \ A\3f „ „-„N \ 



V 



v ud 




v ub 


Vcd 


Vcs 


v cb 


V td 


V ts 





J 



A 2 /2 
-A 



. A\ 3 (l - p-irj) 



A AX 3 (p-ir]) 
1 - A 2 /2 AA 2 
-AA 2 1 



(29) 



for the Wolfenstein parameters appearing in the above expression, we shall use the values given 
by [201 

A = 0.2236, A = 0.824, p = 0.22, fj = 0.35, (30) 

where p = p(l — 4r) and fj = 7/(1 — 4r). For computing the branching ratio, the lifetime of B s 
meson is t Bs = 1.461 ps |2"Uj . 

For the LCDAs of the vector meson, we use the asymptotic form [Fj\ 



(a)V 

91 



6x(l — x) 



4 



[1 + (2x 



(31) 



As for the two B s meson wave functions given by Eq. (jl8|) . we find that only $f (£) has con- 
tributions to the nonfactorizable corrections. We adopt the moments of the $f (£) defined by 
|Hl E] in our numerical evaluation 



(32) 

o 4 A B 

with A# = 0.35 GeV. The quantity A# parameterizes our ignorance about the B s meson 
distribution amplitudes and thus brings large theoretical uncertainty. 

The decay constants and form factors are nonperturbative parameters which are taken as 
input parameters. In principle, they are available from the experimental data and /or estimated 
with well-founded theories, such as lattice calculations, QCD sum rules etc. For the decay 
constants, we take their values in our calculations as El EU E] 



f Bs = 236 MeV, f K * = 214 MeV, f£* = 175 MeV, f„ = 210 MeV, 



195 MeV, 



233 MeV, f± = 175 MeV. 



(33) 



For the form factors involving the B s — > K* and B s — > <p transition, we adopt the results given 
by ^1] which are analyzed using the light-cone sum rule(LCSR) method with the parameteri- 
zation 

/(°) 



1 - a F (q 2 /m Bs ) + b F (q 2 /m Bs ) 2 
10 



(34) 



for the form-factor q 2 dependence. At the maximum recoil, the form factors are listed as ^3] 



A B s<t>(C\\ 


— u.zyo, 


dp = 


U.o 1 , 


Op — 


— U.Uol, 


A 2 l U J 


— u.zoo, 


cip = 


l.OO, 


u 

Op — 


U.olo, 


V B °<t>(o) 


= 0.433, 


of = 


1.75, 


bp = 


0.736, 


Af° K *(0) 


= 0.190, 


dp = 


1.02, 


bp = 


-0.037, 




= 0.164, 


ap = 


1.77, 


bp = 


0.729, 




= 0.262, 


dp = 


1.89, 


bp = 


0.846. 



It is obvious that the q 2 dependence for the form factors A 2 and V are dominated by the dipole 
terms, while A\ by the monopole term in the region where q 2 is not too large. 

To illustrate the non-universality of the nonfactorizable effects on different helicity ampli- 
tudes, we list a few numerical results of the parameters for a specific mode B s K + *p~ in 
Table 1. In order to compare with the parameters di in the NF approach, we also present the 
results of di calculated in NF approach. 

Table 1: The effective parameters d\ in the NF and QCDF approach for B s — > K + * p~ . 





NF 


QCDF 


d\ 


1.0193 


1.0265 + 0.0126i 


a\ 


-0.0293 


-0.0263 - 0.0015i 


d° 

u 10 


-0.0013 


-0.0009 + 0.0007i 


a f 


1.0193 


1.0701 + 0.0126i 


dt 


-0.0293 


-0.0385 - 0.0015* 




-0.0013 


0.0015 + 0.0007z 


d^ 


1.0193 


1.0943 + 0.0126z 


a 7 


-0.0293 


-0.0374 + 0.0022i 


a w 


-0.0013 


0.0028 + 0.0007z 



From Table 1, we can see that nonfactorizable corrections to the helicity amplitudes are 
not universal. The effective parameters d\ for helicity h = 0, +, — states receive different 
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nonfactorizable contributions and hence they are helicity dependent, quite contrary to the naive 
factorization(NF) approach where the parameters are universal and polarization independent. 

The branching ratios for several channels of B s — > VV decays in the LCSR analysis for form 
factors are collected in Table 2. In order to compare the size of different helicity amplitudes, 
we define two quantities: 

rv 
r 



r 







| 2 + 


\ H - 


12 




2 + 


\H+ 
| 2 + 


| 2 + 
\H_ 


\H-\ 2 
12 


\H 


2 + 


\H + 


l 2 + 


\ H -\ 2 



(36) 
(37) 



The ratios of IV/r and T^jT measure the relative amount of transversely and longitudinally 
polarized vector meson. In Table II, we also give the values of IV/r for each channel both in 
QCD factorization(QCDF) approach and the naive factorization(NF) approach. 

Table 2: Branching ratios and the transverse to total decay rate IV/r for charmless B s — > VV 
decays in QCD factorization(QCDF) approach and in the NF approach. 



channel 



r T /r BR 



QCDF NF QCDF NF 



b — > d transition 

B s -> K + *p- 0.071 0.066 1.82 x 10~ 5 1.79 x 10~ 5 

B s -> K°*p° 0.072 0.062 5.29 x 10" 7 5.94 x 10~ 7 

B s -> K°*u 0.046 0.064 3.08 x 10~ 7 7.32 x 10~ 7 
b — > s transition 

B s -> K + *K~* 0.100 0.103 1.94 x 10" 6 1.58 x 10~ 6 

B s -> u!(p 0.141 0.072 5.31 x 10~ 7 2.64 x 10~ 7 

B s -> p°0 0.089 0.070 1.03 x 10~ 6 7.14 x 10~ 7 



pure penguin processes 

B s -> K°*K°* 0.094 0.082 2.61 x 10~ 6 1.94 x 10~ 6 

B s ^K°*(f) 0.190 0.092 1.35 x 10~ 7 1.41 x 10~ 7 

B,^d>d> 0.134 0.117 1.31 x 10~ 5 9.05 x 10~ 6 



From Table 2, we can find that some channels have large branching ratios of order 10 
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which are measurable at near future experiments at CERN LHCb. Owing to the absence of 
(S — P)(S + P) penguin operator contributions to iy-emission amplitudes, tree-dominated 
B s — > V1V2 decays tend to have larger branching ratios than the penguin-dominated ones. 
Moreover, we find that the transverse to total decay rate T^/r is very small for most decay 
modes, so in the heavy quark limit, both light vector mesons in charmless B s —>■ VV decays 
tend to have zero helicity. 

4 Summary 

In this paper, we calculated the branching ratios for two-body charmless hadronic B s —* VV 
decays within the framework of QCD factorization. Contrary to phenomenological generalized 
factorization [9] and NF approach, the nonfactorizable corrections to each helicity amplitude 
are not the same. The effective parameters at? vary for different helicity amplitude and hence 
are helicity dependent. Since the leading-twist DAs contributions to the transversely polarized 
amplitudes vanish in the chiral limit, in order to have renormalization scale and scheme inde- 
pendent predictions, it is necessary to take into account the contributions of the twist-3 DAs 
of the vector meson. Contrary to the PP and PV modes, the annihilation amplitudes in the 
VV case do not gain the chiral enhancement of order mt B j im q m^). So we do not include the 
contributions of the annihilation diagrams which is truly power suppressed in the heavy quark 
limit. It should be stressed that we have not taken into account the higher-twist DAs contri- 
bution for the longitudinally polarized vector meson. Though direct calculation, the transverse 
to total decay rate IV/r is found to be very small, and both light vector mesons tend to have 
zero helicity. Branching ratios of B s — > VV decays are calculated with the LCSR analysis for 
the form factors and the branching ratios of some channels are found as large as 10~ 5 , which 
might be accessible at future experiments at CERN LHCb. 
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Appendix. 

The B s — > VV decay amplitudes are collected here: 
I. b — > d processes: 



rhrrs 











■* K°*lu) 




, and Af 


= 



3 h 3 
-aS + -< 
2 7 2 



H\B s ^p^) = ^{X u 4-\ t [ 3 -(a h 7+ a h 9 )]}x^ 
H h (B.^u<f>) = ^^X u a h 2 -X t [2(a h 3 + al) + ^ + a^X^\ 

where X u = V ub V* s and X t = V tb V t * s 
3. pure penguin processes: 



G, 



H h (B s — > K°*(p) = -^ W {[4 + ^-i(4 + 4)]X^°^) 

H h (B s ^W) = -^ Vtb V;Aa h 3 + 4 + a h 5 -^ + a h 9 +a h w )]X^ 
In the above expressions, the factorozable amplitude X^ BaVl ' V2 ^ is defined as in Eq 
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